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Abstract 

The natural orbitals and natural occupation numbers of various N = Z, sp and sd shell nuclei 
are calculated by applying a correlated one-body density matrix. The correlated density matrix has 
been evaluated by considering central correlations of Jastrow type and an approximation named 
factor cluster expansion. The correlation effects on the natural orbitals, natural occupation numbers 
and the Fermi sea depletion are discussed and analysed. In addition, an approximate expression for 
the correlated one-body density matrix of the nuclear matter has been used for the evaluation of the 
relative momentum distribution and the Fermi sea depletion. We found that the value of the Fermi 
sea depletion is higher in closed shell nuclei compared to open shell ones and it is lower compared 
to the case of nuclear matter. This statement could be confirmed by relevant experimental studies. 
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I. INTRODUCTION 



The short-range correlations (SRC) play an important role on the one- and two-body 
properties of nuclei and nuclear matter. In general, the contribution of SRC is important 
for the description of the mean value of some two-body operators, such as the ground- 
state energy of nuclei, but it is also of interest to investigate the SRC contribution to simpler 
nuclear quantities related to one-body operators such as the form factor, density distribution 
and momentum distribution. The key quantity to calculate the one-body properties of a 
quantum many-body system is the one-body density matrix (OBDM)[l|. The knowledge 
of the OBDM leads to the evaluation of the one-body properties of a quantum many-body 
system both in coordinate and in momentum space. In addition, the diagonalization of the 
OBDM leads to the evaluation of the natural orbitals (NO) and natural occupation numbers 
(NON). Over the last decades there was an effort to incorporate SRC on the OBDM mainly 
for light and medium nuclei [2 17]. 

The OBDM has been evaluated recently, by considering central correlations of Jastrow 



type 



181 ] in the framework of the Iwamoto and Yamada factor (FIY) cluster expansion [19] 



developed by Clark and Westhaus [20| and Feenberg |2Jj. More precisely, the expression of 
;he OBDM, p(r, r'), was found by using the factor cluster expansion of Clark and co-workers 
22j and Jastrow correlation function which introduces SRC for N = Z sp and sd shell nuclei. 
The evaluated OBDM is a functional of the harmonic oscillator (HO) orbitals and depends 
on the HO parameter b = (h/mu) 1 ^ 2 and the correlation parameter (3 |S|] . 

In addition, by employing a phenomenological ansantz for the correlated OBDM of the 
uniform nuclear matter we calculated the relative momentum distribution and consequently 
the Fermi sea depletion (FSD). We considered a value for the nuclear matter density close 
to equilibrium density of symmetrical nuclear matter in order to be able to compare the 
correlated parameter and consequently the effect of SRC in finite nuclei and in nuclear 
matter. 

The motivation of the present work is twofold. Firstly, we applied the OBDM correspond- 
ing to finite nuclei in order to make a detailed investigation of the correlation effects on the 
NO, the natural occupation probabilities and the depletion of the Fermi sea in light and 
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medium nuclei. Actually, very little information can be found in physics literature about the 
properties of natural orbitals and occupation probabilities. In view of the above statement, 
we additionally extend the study of the NO and NON not only to closed shell but also to 
some open shell nuclei. To our knowledge this is an intriguing study that has never been 
examined. Secondly, we extended the calculations in nuclear matter in order to compare the 
effect of the SRC in finite and infinite nuclear matter under almost the same conditions and 
to provide results which could pose for future experimental investigation. 

The article is organized as follows. In Sec. II the correlated OBDM and the diagonalization 
process in finite nuclei are presented, while the relative theory of nuclear matter is reported 
in Sec. III. The results are discussed in Sec. IV. Finally the summary of the work is given in 
Sec. V. 

II. CORRELATED ONE-BODY DENSITY MATRIX AND NATURAL ORBITALS 
REPRESENTATION 

A. One-body density matrix 

A nucleus with A nucleons is described by the wave function ^(rj, r 2 , ■ • • , r^) which 
depends on 3A coordinates as well as on spins and isospins. The evaluation of the single 
particle characteristic of the system requires the one-body density matrix [l| 



where the integration is carried out over the radius vectors r 2 , • • • , and summation over 
spin and isospin variables is implied. p(r, r') can also be represented by the form 



where ty' = r' 2 , ■ ■ • , r^) and iV is the normalization factor. The one-body density 

operator O rr >(A), has the form 




(1) 



p(r, r') 



N(y\O rr ,{A)\V) = N(O rr ,{A)), 



(2) 



.4 



A 




(3) 
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The diagonal elements of the OBDM give the density distribution p(r, r) = p(r), while 
the momentum distribution is given by the Fourier transform of p(r, r'), 

nCk) = -. — rr / exp[zk(r — r')]p(r, r')drdr'. (4) 
(2vr) d J 

The trial wave function which describes a correlated nuclear system, can be written as 

* = J 7 ®, (5) 

where $ is a model wave function that is adequate to describe the uncorrelated A-particle 
nuclear system and J 7 is the operator that introduces SRC. The function $, is chosen to be 
a Slater determinant wave function, constructed by A ortho-normalised single-particle wave 
functions </>i(r) 

$ = $ SD (n, r 2 , • ■ ■ , v A ) = (v4!)~ 1/2 det | ^(r,) |, (i, j = 1, 2, • • • , A). (6) 

Several restrictions can be applied on the model operator J 7 . In the present work T is taken 
to be of Jastrow type [18 1 

A 

i<j 

where /(ry) is the state- independent correlation function of the form 

/(r iJ ) = l-exp[-/3(r i -r J ) 2 ]. (8) 

The correlation function /(r^-) tends to 1 for large values of = |r< — | and it tends to 
for — ► 1. 

The expression of the correlated one-body density matrix p(r, r') has been found by 
applying the method of the generalised integrals and by using the factor cluster expansion 
of Ristig, Ter Low, and Clark j^]. The various terms of p(r, r') are given in Ref. js]. 
The relative expression depends on the single-particle wave functions and so it is suitable 
to be used for analytical calculations with the HO orbitals and in principle for numerical 
calculations with more realistic single-particle orbitals. The calculations were carried out for 
sp and sd closed shell and open shell nuclei with N = Z. 
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B. Mean-field approximation and harmonic oscillator wave functions 

In the mean field approximation (MFA) the single-particle wave functions 4> a (r) appearing 
in the A-body wave function, given by Eq. ([6]), form a nuclear Fermi sea F. These wave 
functions are called hole-state orbitals and define the uncorrelated OBDM 

^^smw- ( 9 ) 

Eq. fl9]) is associated with the Slater determinant wave function given by Eq. (|6]). The 
uncorrelated density matrix pmfa{ y , v ') is diagonal and the occupation probabilities are 
equal to unity inside the Fermi sea and zero for the states (called particle states) outside 
the Fermi sea. In the present work the Pmfa{ v i r ') was constructed by using single-particle 
wave functions that originated from the HO potential. 



C. Natural orbitals representation of the OBDM 

A model-independent way to define a set of single-particle wave functions and occupation 
probabilities from the correlated p(r, r') is to use the natural orbital representation suggested, 
by Lowdin [1], according to the following expansion 

p(r,r') = 5> a C(r)^(r'). (10) 

a 

The normalized eigenfunctions ip a (r) °f p( r ) r ') are called natural orbitals and they form 
a complete orthogonal set. The eigenvalues n a , where < n a < 1 are called natural 
occupation numbers and they define the occupation probability of the NO ip a {r). Usually 
there are orbitals i/j a (r) for which the occupation numbers n a are significantly larger than 
the others. These are called hole state orbitals and form a new Fermi sea, while the rest are 
called particle state orbitals Q]. 

The density distribution in natural orbitals representation is of the form 

p(r) = 5> Q |^(r)| 2 . (11) 
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The normalization is satisfied via the sum rule 



^n a = A or y^p Q = 1, (12) 

a a 

where p a = n a /A is the occupation probability of the state a. 

We apply the one-body density matrix, evaluated for the various sp and sd shell nuclei js]. 
The various terms of p(r, r') are functionals of the harmonic oscillator orbitals depending on 
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the size parameter b and the correlation parameter (or y = fib 2 ) 

In general, the diagonalization of p(r,r'), according to definition fflO|) . is achieved by 
solving the equation 

J p(r,v')ip a (r')dr' = p a ?p a (r). (13) 

For nuclei with total angular momentum J = 0, the OBDM has to be diagonalized with 
the {Ijm} subspace of the complete space of NO. For the case of the closed shell nuclei their 
angular part is explicitly determined 

VUr) = Vw(r) = iJ nl (r)Yr($, <p). (14) 

Substituting Eq. (fT4|) into f|T3|) and integrating over the angular variables, the equation for 
the radial part of the NO becomes 

47T f°° 

2l - ~ J pi(r,r')tp nl (r r y 2 dr' = p nl tp nl (r). (15) 
The kernels pi(r,r') are the coefficients of the Legendre expansion of the OBDM 

oo 

p{r,r') = p(r,r', cos u rr >) = 2^Pi(r, r')Pi (cos u rr i), (16) 

1=0 

where 

21 + 1 Z" 1 

pi(r,r') = / p(r,r' , cos u rr /)Pi(coscu rr r)d(coscu rr >). (17) 

2 J-i 

From the analytical expression of the OBDM, p(r, r') = p(r, r', cos u) rr '), the analytical 
expression of the expansion coefficients pi(r,r') were found and they have the form [8| 

3 

pi (r, r') ~ iV [p u (r, r') - p 2 2« (r, r', ^ 

i=l 
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where the one-body contribution of the expansion coefficients and the three terms 
P22i{r, r', gi) (i = 1,2,3), which come from the two-body contribution of the cluster ex- 
pansion, have the form 

pu(r,r') = * [(4r] ls +r] 2s U Q (r b , r' b ))6i +rii p n. 1 (r b , r b )5n+r] ld U 2 (r b , r b )Si 2 ] exp[-(rl+r' 2 b )/2], 

(19) 

P22/ (r, r , gi )= exp [- — - - r b - r b /2] 

x [A i(r b ,r' b ,gi)5io + Au(r b ,r' b , g^Sn + A 2 i(r b ,r' b , &)($o + 25 i2 )/3], % = 1,2 

(20) 

2/ + 1 4 

P22/(r, r', c/ 3 ) = — ^— B w (r 6 , r£, g 3 ) [fuiin, r' b )e ci + f 1M (r b , r' b )e C2 ] . (21) 

fc=0 

In Eqs. ([19]), (J2DJ) and ((21]) n i; / 2W , ci and c 2 are polynomials of r b = r/b 

and r' b = r'/b. It is noted that pu(r,r') = P22i( r , r ' , 9i) = P22i( r , r ' , 92) = for / > 2, while 
P22i(r,r',g 3 ) ^ V/. 

Using the analytical expressions of pi(r,r') we solved the eigenvalue equation ( fl5l) and 
as a result the dependence of the occupation probabilities p n i and the NO ipniij') on the 
parameters b and y can be studied. 



III. ONE BODY DENSITY MATRIX AND MOMENTUM DISTRIBUTION IN 
NUCLEAR MATTER 

The model we used is based on the Jastrow ansatz for the ground state wave function of 
nuclear matter 

#(ri,ra,...,rjv)= JJ f(r ij )^ (v 1 ,r 2 ,...,r N ), (22) 

l<i<j<N 

where = |r$ — | , $0 is a Slater determinant (here constructed by plane waves with 
appropriate spin-isospin factors, filling the Fermi sea) and f(r) is a state- independent two- 
body correlation function. 

A cluster expansion for the one-body density matrix p(r, r') of the nuclear matter, in the 
framework of the low order approximation (LOA), has been derived by Gaudin, Gillespie 
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and Ripka [25|, |26] for the Jastrow trial function (1221) and has the form 

p LOA (r,r') = p (r,r') + J [f(\r - r 2 |)/(|r' - r 2 |) - 1] 

x [4p (r 2 , r 2 )p (r, r') - p (r, r 2 )p (r 2 , r')] dr 2 

~ J J [/ 2 (|r 2 -r 3 |)-l]po(r,r 2 ) 

x [4p (r 2 , r')p (r 3 , r 3 ) - p (r 2 , r 3 )p (r 3 , r')] dr 2 dr 3 , 



(23) 



where po(r, r') is the density matrix corresponding to the wave function $o (that is f(r) = 1). 
Formula (1231) is also applicable to finite nuclei j^, 6 1 . In the case of uniform and infinite nuclear 



matter p 



that 
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_ 1 QA(r,r / ) becomes only a function of u = |r — r'|. In addition it is easy to prove 
26] 



Ploa(m) 
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(2vr) : 



e~ iku dk. 



(24) 



Obviously Ploa(O) = Po, where p is the density of the nuclear matter. The corresponding 
momentum distribution is given by the Fourier transform 

1 



^loaO) = \ J PhOA(u)e tku du. 
Adopting the gaussian model for the correlation function f(r) 

/(r) = l-exp[-/3r 2 ], 



the momentum distribution can be written 26] 

n hQ A(k) = 9{kp — k)[l — k diT + Y(k, 8)] + 8 [k dir Y(k, 2) - [Y(k, 4)] 2 ] 



where 



and 



c- l Y{k,p) 



-P -k 2 

e + — e - 
2k 



2 ~ 

e~ y dy + sgn(A;_ 



|fe_ 



dy, 



1 



p\3/2 



k 



k 



k, 



k F ± k 



p = 2,4,8, 



(25) 
(26) 

(27) 

(28) 
(29) 



while sgn(x) = 

The dimensionless Jastrow wound parameter fcdir, which can serve as a rough measure of 
correlations and the rate of convergence of the cluster expansion, is defined as 



^dir = Po J [f(r) - l] 2 dr, 



(30) 



where po is the density of nuclear matter, which is related to the Fermi wave number fcp by 
the equation po — 2/cf 7 /(37r 2 ). The normalization condition for the momentum distribution 
is 

I 

n LOA (k)k 2 dk = - k F . (31) 

From Eq. ( 1301) we obtain the following relation between the wound parameter k^h and the 
correlation parameter 

1 / k, x 4 



k - = *m\w- (32) 

It is clear that large values of fcdir imply strong correlations and poor convergence of 
the cluster expansion. A reasonable interval of the values of the correlation parameter 
(3 is: 1.13 fm~ 2 < < 7 fnT 2 ^26\. That range corresponds to 0.3 > k^h > 0.02. In our 
calculations we considered the value po = 0.182 fm -3 for the density of uniform, spin-isospin- 
saturated nuclear matter which corresponds to kp = 1.3915 fm -1 . This value is close to hp 
of symmetrical nuclear matter at equilibrium density. It is mentioned that Eq. (T23]) and 
(I27|) which are applicable to nuclear matter at lower effective densities, can also be applied 



reliably to nuclei 
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26]. Thus, a comparison of the correlated parameter i.e. of the effect of 



SRC in finite nuclei and in nuclear matter can be made by considering a value of the nuclear 
matter density close to equilibrium density of symmetrical nuclear matter. 



IV. RESULTS AND DISCUSSION 

I — I 

We used the one-body density matrix, calculated in our previous work [8j, in order to 
construct a corresponding set of natural orbitals and natural occupation numbers for various 
sp and sd shell nuclei, with N = Z. In that model there are two free parameters, the size 
parameter b and the correlation parameter /3 (or y = j3b 2 ). These have been determined by 
least squares fit to the experimental form factor and are given in Table I of Ref. 8j. It is 
worth noting that both the NO and the NON depend strongly on the correlation parameter 
y and one may study extensively the effect of the SRC on the above quantities. However 
the purpose of the present work is to use fixed parameters from the previous work and to 
compare NO and NON with other similar studies and the existing experimental data. 
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First we define the occupation ratio r\ a of the state a as 

number of nucleons occupy the state a . . 

maximum capacity of the state a 

The ratio r\ a is related to the occupation probability p a according to 

_ Pa ' A 

01 maximum capacity of the state a 
In view of the above definition the occupation ratios of the various NO states will be given 
by the expressions 

Ps-A p p -A p d - A p r A pg-A 

Vs = —, Vp = —, m = — ^ = -gg- (35) 

The values of the occupation ratio r\ a of various states for various sp and sd shell nuclei are 
presented in Table [B It is noted that in the MFA (absence of correlations) the corresponding 
value 7/^ FA for the closed shell nuclei is 1 for the hole states and for the particle states. 
For non closed shell nuclei, we considered that for the MFA hole states lp for 12 C and Id 
for 24 Mg, 28 Si, 32 S and 36 Ar the occupation ratios r/^ FA are: 0.6667, 0.4, 0.6, 0.8 and 1, 
respectively. 



The values of r] a , given in Table [J are compared with experimental data 27] and the 



It should be noted that the values of the 



theoretical estimate obtained for 40 Ca in Ref. [6]. 
occupation ratios r] a are sensitive on the various elements of the correlated OBDM as well 
as on the strength of the correlations. The use of different methods in the evaluation of the 
OBDM may lead to quantitatively different values of 7] a . 

In a few cases arose non-physical negative values for some occupation ratios r] a while 
in two cases (Id and 2s states of 40 Ca) the occupation ratios rjid and r]2 S surpass slightly 
the upper limit, that is r)id,T)2s > 1- Those results originated from the breaking of the 



/L-representability [28J of the OBDM applied in the present work. The breaking of the A- 
representability depends both on the form of the correlated one-body density matrix as well 
as on the strength of the correlations which were introduced via the correlation parameter 
y. In general the breaking of the A- representability expresses, in a quantitative as well as 
qualitative way, the pathology of the approximations made in the expansion of the density 
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matrix p(r, r'). One may be able to restore this problem by excluding all states with negative 
f] a and eventually renormalize the resulting A-representable density matrix p(r, r') [6] but 
this procedure is out of the purpose of the present work. 
The depletion of the hole states, which is defined as 

D a = (l - Jp^) • 100 (%) , (36) 
and the Fermi sea depletion Dp, which is defined as 

Df = j (VTA* + 12C A ^i P + 20tjff A D u + ^ FA D 2s ) (%) , (37) 
are given in Table [TT1 

It is instructive to compare the present values of the FSD with those originating from 
other theoretical and experimental studies. Actually, the FSD, compared with the specific 
depletion of the various states, can be measured, at least indirectly or with model depen- 
dent analysis, in relative experiments. The value of the nuclear FSD, in the present work, 
according to Table [His: 5.16% for 4 He, 3.41% for 12 C, 4.00% for 16 0, 1.59% for 24 Mg, 2.85% 
for 28 Si, 2.52% for 32 S, 5.38% for 36 Ar and 5.91% for 40 Ca. It is obvious that the value of 
the FSD is higher in the closed shell nuclei 4 He, 16 and 40 Ca where the SRC are stronger 
compared to the neighbourly open shell nuclei. The above statement can only be confirmed 
by future experimental studies of the Fermi sea depletion which will include both closed and 
open shell nuclei. The present results concerning the FSD of closed shell nuclei, are com- 



parable to those evaluated by Stoitsov et al. [6J. In the mentioned reference the depletions 
of the Fermi sea were 5.7%, 3.9%, and 6.3% for 4 He, 16 and 40 Ca, respectively. However, 
for the case of 40 Ca, the depletion is smaller compared to the value 9.4% obtained in the 
analysis of the experimental data from (e, e'p) and (d, 3 He) reactions on the 40 Ca nucleus 



27] 



In Figure Q] and [2} we compare the uncorrelated HO single-particle wave functions with 
the corresponding correlated natural orbitals for the considered nuclei. The most distinctive 
feature is the similarity of the HO and the NO single-particle wave functions for the states 
lying below Fermi level. This is an indication that the effect of short-range correlations on 
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the configuration space orbitals is not significant for nucleon moving below Fermi level. On 
the contrary, the main difference between the HO and the NO wave functions appears for 
the states lying above the Fermi level. Due to SRC, the particle NO wave functions are 
concentrated mostly inside the nucleus while the HO ones are expanded outside the nucleus. 
As a result the natural particle-state orbitals have smaller rms radii and significantly larger 
high-momentum components compared to the unoccupied MFA orbitals. In addition, they 
have small but non-zero occupation ratios. The above properties of the natural orbitals 
affects appreciably the nucleon momentum distribution n(k). More precisely it was found 
that the high momentum components of the momentum distribution n(k) can be determined 
to a large extent by the particle-state contribution due to SRC (| . 

The behaviour of the momentum distribution of the uniform nuclear matter, by employing 
Eq. (j27p . as a function of kjkp for various values of the wound parameter k^ r is indicated in 
Figure [3](a). The discontinuity Zp, of the momentum distribution at kjkp = 1, is another 
characteristic quantity used as a measure of the strength of correlations of the uniform Fermi 
systems and is defined as 

Z F = n(l~) -n(l + ). (38) 

The discontinuity Zp for various values of fcdir is displayed in Figure [3]^a). Zp = 1 for 
the uncorrelated matter, while for correlated matter Zp < 1. In the limit of very strong 
correlations Zp = there is no discontinuity on the momentum distribution of the nuclear 
matter. The quantity (1 — Zp) measures the ability of correlations to deplete the Fermi sea 
by exciting particles from the hole states to the particle states 26]. 

The Fermi sea depletion D nm in the case of uniform nuclear matter is defined as 



D nm = (^3 J n LOA (x)x 2 dxj ■ 100 (%), x = kjkp. (39) 

The depletion of the Fermi sea of finite nuclei and nuclear matter as a function of the 
correlation parameter is shown in Figure [3]^b). As a comparison we also include the results 
of Ref. 6j where the OBDM of finite nuclei was evaluated using Eq. ( [231) . This equation was 
employed in the present work for the study of the Fermi sea depletion of nuclear matter. In 
a similar work, Jaminon et al. used the relativistic Brueckner-Hartree-Fock approximation 
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and they found that the FSD is around 5% [29]. In addition, de Jong et al, in the framework 
of a similar relativistic model found that the average depletion of the Fermi sea is 11% [30j. 

It is obvious that the value of the FSD, for the same value of the correlation parameter (3 is 
higher in nuclear matter compared to finite nuclei. Actually, even for the same approximation 
for the OBDM (that is the LOA given by Eq. (J23l) ) the depletion in nuclear matter is almost 
twice compared to finite nuclei. However it should be noted that the dependence of FSD on 
j3 exhibits almost similar behavior for nuclear matter and finite nuclei. This is an indication 
that although the finite nuclei have a different structure compared to uniform nuclear matter, 
SRC affect their mean field structure (uncorrelated case) in the same manner at least for 
nuclear matter density close to the saturation density. 



V. SUMMARY 



We studied the effect of short range correlations on the natural orbitals and natural occu- 
pation numbers for various sp and sd shell nuclei via the correlated one-body density matrix. 
A comparison with the MFA is also presented within the framework of the Fermi sea de- 
pletion. It was found that the effect of the SRC on the natural orbitals is not significant 
for nucleons moving below the Fermi level. The main difference between HO and NO wave 
functions appears for states lying above the Fermi level affecting appreciably the high mo- 
mentum components of the momentum distribution n(k). Finally, we found that the value 
of the Fermi sea depletion is higher in the closed shell nuclei compared to the open shell 
nuclei but in both cases it is lower compared to that of nuclear matter even for the same 
strength of SRC. However definitive conclusions can only be obtained by future experimental 
investigation. 
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FIG. 1: The natural orbitals of the particle- and hole-states (solid lines) compared to the cor- 
responding HO orbitals (dashed lines) for the closed shell nuclei 4 He, 16 0, 40 Ca as well as for 
36 Ar. 
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FIG. 2: The natural orbitals of the particle- and hole-states (solid lines) compared to the corre- 
sponding HO orbitals (dashed lines) for the open shell nuclei 12 C, 24 Mg, 28 Si and 32 S. 
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FIG. 3: (a) The momentum distribution for correlated nuclear matter versus kjkp for various 
values of the correlation parameter (the values of n(x) for x > 1 have been multiplied by 10) 
(b) The Fermi sea depletion of nuclear matter (solid line), and finite nuclei of the present work 
(solid circles) and of Ref. [g] (open circles). 
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TABLE I: The occupation ratio rj a of particle- and hole-states for various sp and sd shell nuclei. The 
last two columns correspond to theoretical and experimental data of Ref. [g] and |27J| , respectively. 



nl 4 He 12 C 16 24 Mg 28 Si 32 S 36 Ar 40 Ca 40 Ca [6j 40 Ca [27] 

Is 0.9484 0.9044 0.8846 0.9497 0.9066 0.9200 0.8156 0.6630 0.8899 0.9500 

2s -0.0037 0.0029 0.0036 0.0009 0.0023 0.0019 0.0066 1.0686 0.9588 0.8900 

3s 0.0037 -0.0017 -0.0023 0.0005 0.0012 0.0011 0.0030 0.0089 -0.0097 

4s 0.0002 0.0009 0.0013 0.0003 0.0005 0.0005 -0.0017 -0.0089 0.0113 

lp 0.0077 0.6644 0.9852 0.9794 0.9564 0.9580 0.9021 0.8889 0.9296 0.9200 

2p 0.0008 0.0026 0.0031 0.0008 0.0020 0.0016 0.0053 0.0093 -0.0055 0.1200 

3p 0.0001 0.0007 -0.0022 0.0004 0.0010 0.0009 -0.0039 -0.0054 0.0121 

4p 0.0000 -0.0005 0.0009 0.0002 -0.0008 -0.0005 0.0022 0.0037 0.0041 

Id 0.0025 0.0036 0.0052 0.4034 0.6050 0.8009 0.9988 1.0022 0.9468 0.8920 

2d 0.0003 0.0011 0.0017 0.0008 0.0020 0.0016 0.0047 0.0066 -0.0063 

3d 0.0000 0.0003 0.0005 0.0004 0.0009 0.0008 0.0042 -0.0052 0.0082 

4d 0.0000 0.0000 0.0002 0.0002 0.0003 0.0003 0.0018 0.0027 0.0028 

If 0.0008 0.0020 0.0033 0.0009 0.0025 0.0024 0.0080 0.0090 0.0127 0.3700 

2f 0.0001 0.0006 0.0010 0.0005 0.0011 0.0010 0.0031 0.0041 0.0050 

3f 0.0000 0.0002 0.0003 0.0002 0.0005 0.0005 0.0012 0.0017 0.0018 

4f 0.0000 0.0000 0.0001 0.0001 0.0002 0.0002 0.0004 0.0007 0.0006 

lg 0.0003 0.0013 0.0021 0.0007 0.0019 0.0017 0.0059 0.0065 0.0087 

2g 0.0000 0.0004 0.0006 0.0003 0.0008 0.0008 0.0022 0.0028 0.0031 

3g 0.0000 0.0001 0.0002 0.0002 0.0003 0.0003 0.0008 0.0011 0.0010 

4g 0.0000 0.0000 0.0000 0.0001 0.0001 0.0002 0.0003 0.0004 
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TABLE II: The depletion of the hole states and the FSD (in %) for various sp and sd nuclei 
compared to the theoretical and experimental data of Ref . |6fl and 27[ , respectively. 



nl 4 He 12 C 16 24 Mg 28 Si 32 S 36 Ar 40 Ca 40 Ca [6] 40 Ca [27] 

Is 5.16 9.56 11.54 5.03 9.34 8.02 18.44 33.7 11.01 5.00 
2s -6.86 4.12 11.00 

lp 0.34 1.48 2.06 4.36 4.20 9.79 11.11 7.04 8.00 

Id -0.85 -0.83 -0.11 0.12 -0.22 5.32 10.80 

FSD 5.16 3.41 4.00 1.59 2.85 2.52 5.38 5.91 6.30 9.40 
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